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ABSTRACT 


The Apollo Lunar Program utilized efficient, t.c., Earth-return, transearth 
trajectories which employed parking orbits in order to minimize energy requirements. 
This thesis concentrates on a different type of transearth trajectory. These are direct- 
ascent, hyperbolic trajectories which omit the parking orbits in order to achieve short 
flight times to and from a future lunar base. The object of the thesis is the development 
of a three-dimensional transearth trajectory model and associated computer program for 
exploring trade-offs between llight-time and energy, given various misslon constraints. 
The program also targets the Moon with a hyperbolic trajectory, which can with a time- 
reversed trajectory; be used for targeting Earth impact points. The first-order model is 
based on an Earth-centered conic and a massless spherical Moon, using MathCAD 
version 3.0. This model is intended as the basis for future patched-conic formulation for 
the design of fast Earth-return trajectories. Applications tnclude placing nuclear- 
deterrent arsenals on the Moon, various space support related activities and finally 


protection against Earth-threatening asteroids and comets using lunar bases. 
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TABLE OF SYMBOLS AND ABBREVIATIONS 


This listing of symbols and abbreviations is to aid the reader in clarifying the 


multitude of symbols and abbreviations used in this thesis. They are in alphabetical order. 


Or ORB 


Ср 


An? 
QEM 
Сұм] 
QEM? 


QEN 


Со 


Ве 


: input value for posigrade or retrograde orbit 
: spacecraft right ascension at Earth launch 

: spacecraft right ascension at Moon intercept 
: represents cosine of a, 

: represents sine of a, 

: represents the difference between a, and ag 
: represents the cosine of agy 

: represents the sinc olan 


: represents the difference between the translunar trajectory longitude-of- 


ascending-node and spacecraft right ascension at Earth launch 


‚represents the cosine Ol arn 
: represents the sinc Olar, 


: represents an intermediate value to calculate the spacecraft right 


ascension on its trajectory 


: represents spacecraft right ascension on its trajectory 
: represents the selenographic longitude for lunar surface impact 


: represents flight path azimuth at launch 


vi 


Pr 
Br 
Әм 

Pmi 
Pw 


Үні 
УЕ? 
Ум 
У мі 
У м2 
Yo 


Yı 


represents 
represents 
: represents 
горео 5 
: represents 
: represents 
> represents 
; represents 
: represents 
; represents 
; represents 
; represents 
: represents 
; Tepresents 
: represents 
: represents 
тергезеп 
> represents 
> represents 
: represents 
> represents 
‚ represents 


: represents 


cosine of f, 

sine of B, 

[light path azimuth at moon intercept 

cosine of B,, 

sine of By, 

flight path azimuth of the spacecraft on its trajectory 
cosine of Ba 

sine of f, 

the sweep angle at launch 

the cosine of C, 

[Ве ше о С 

the sweep angle at moon intercept 

the sweep angle of the spacecraft on its trajectory 
flight path angle at launch 

the cosine of Yg 

the sine of yp 

light path angle at moon intercept 

the cosine of y,, 

che sine of y 

the flight path angle of the spacecraft on its trajectory 
the cosine of y, 

the sine of y, 


the deflection angle 


Эмо 


Аа 


Дам 


Дам, 
Дау 
Aa, 
Aa, 
Aa, 


А о 


eccen 


COUCHES 


Cnergyai, 


: represents the directrix-to-focus or axial offset of the asymptotic from 


the focus ; it is a variable valid for eccentricities greater than one 


: represents the distance between the spacecraft and the lunar center at 


UME 


: represents the spacecraft declination at earth launch point 

: represents spacecraft declination at moon intercept 

: represents the spacecraft declination on its trajectory 

: represents the selenographic latitude for surface impact on the moon 
«see Figure S 


: represents intermediate value to calculate spacecraft right ascension on 


its trajectory 


: represents the cosine of Aay 

: represents the sinc ol Aaa 

: represents an intermediate value for the calculation of a, 
: represents the cosine oraa; 

: represents the sine of Aa, 

: represents cosine of A, : see iterative process on page 12 


: represents the difference between point-moon intercept flight time 


relative to perigee passage and the spacecraft flight time on its 


trajectory 


: represents eccentricity as an input value 
: represents the minimum eccentricity input value allowed 


: represents the minimum specific energy required for launch 


viii 


: represents the hyperbolic eccentric anomaly for earth launch 
: represents the hyperbolic eccentric anomaly for point-moon intercept 


: represents the hyperbolic eccentric anomaly for the spacecraft on its 


trajectory 


: represents the hyperbolic sine of E,, 
: represents the hyperbolic sine of E,, 
: represents the hyperbolic sine of E,, 


: represents the input value for the true anomaly of the spacecraft at 


moon intercept 


: represents the true anomaly of the spacecraft at earth launch 


: represents the difference of the true anomaly of the spacecraft at moon 


intercept and at earth launch 


: represents a constraint value on fy if required to be used 

: represents the cosine ol fy max 

: represents the sinc of fy мах 

: represents the cosine of f, 

тергезеп the, sine ol f,, 

: represents the truc anomaly of the spacecraft on its trajectory 
: represents the cosine of f, 

: represents the first Euler angle 

: represents the geocentric altitude of launch point 

: represents the inclination of the moon’s orbit plane 


: represents the translunar trajectory inclination 


ІХ 


= 


эл я 


23 
< 


К 


X 
PEM 


518Пєм 


Stem 


: represents 
: represents 
: represents 
: represents 
: represents 
: represents 
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when calculating a, 
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the point-moon intercept flight time relative to perigee 


an intermediate value used to calculate tm 

the point-moon translunar trajectory flight time 
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the difference between the lunar sweep angle at moon 


Intercept and the moon’s mean orbital rotation rate multiplied by A, 


represemi: 


the first Euler angle 


: represents the X-coordinate for the rectangular coordinates of the moon 


at point-moon intercept 


: represents the X-component for calculating radial velocity at point- moon 


intercept 


: represents the X-coordinate for the rectangular coordinates of the 


spacecraft on its trajectory 


: represents the X-coordinate for the center of the moon 
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: represents the Y-coordinate for the rectangular coordinates of the moon 


at point-moon intercept 


: represents the Y-component for calculating radial velocity at point-moon 


intercept 


: represents the Y-coordinate for the rectangular coordinates of the 


spacecraft on its trajectory 


: represents the Y-coordinate for the center of the moon 
: represents the difference between Y, and Y no 


: represents the Y-coordinate for the rectangular selenographic 


coordinates 


: represents the third Euler angle 


: represents the Z-coordinate for the rectangular coordinates of the moon 


at point-moon intercept 


: represents the Z-component for calculating radial velocity at point-moon 
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: represents the Z-coordinate for the center of the moon 
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: represents the Z-coordinate for the rectangular selenographic 
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: represents the Earth's angular rotational rate 
: represents the Moon's mean orbital rotation rate 
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I. INTRODUCTION 


This chapter provides a general description of this thesis and is divided into two 
sections. The first section presents background information pertaining to translunar 
trajectory design. The second section delineates the objectives and limitations of this 


project. 


А. BACKGROUND 

While the impending mission of establishing a permanent lunar base stands in our 
future, it is appropriate to study this type of mission by investigating various translunar 
trajectories. The trade-off between flight time and energy required must be determined 
prior to considering the type of translunar trajectory to design. The energy required for 
the hyperbolic trajectory can then be compared to that of the Apollo mission to 
determine the increase in energy and thereby cost increases. 

The Apollo Lunar Program of the 1960s utilized energy-efficient translunar 
trajectories. This thesis will analyze a different class of trajectories named fast Earth- 
return trajectories. These trajectories utilize a direct-ascent from the lunar surface and 
omit the lunar parking orbit used in the Apollo missions. These direct-ascent trajectories 
differ considerably from translunar trajectories previously considered and discussed in this 


thesis. (Wadsworth, 1991) 


В. OBJECTIVES AND LIMITATIONS 

The objective of this thesis is to develop a trajectory model and associated 
computer program for systems studies of fast Earth-return trajectories from a lunar base. 
The scope of the project involves a step process which evolves in complexity, as the 
design moves from a massless point-Moon, to a massless spherical Moon intercept 
solution. Finally, the groundwork is laid for a patched-conic model which culminates in 
a solution for intercept on a Moon with finite lunar mass. A novel feature is the use of 
a time-reversed trajectory for targeting Earth impact points. In other words, the lunar 
impact point serves as the lunar launch point tor a time-reversed trajectory. The 
complete patched conic approximation, as discussed in Chapter II, is not developed in the 


program due to scope of work and time constraints. 


H. ANALYSIS APPROACII 


This chapter sets the foundation of the thesis by delineating a flow chart of tasks 
that are required to be accomplished to reach the objective of the thesis. A flow chart 
represented in Figure 1, summarizes the work to be accomplished. The last block of the 


diagram represents further rescarch work. 


DESIGN A LUNAR TRAJECTORY FIND THE INTERSECTION OF 
TO INTERCEPT A MASSLESS THE EARTH-CENTERED 
POINT-MOON MOVING IN A CONIC TRAJECTORY WITH 


THE SURFACE OF A 
MASSLESS SPHERICAL MOON 


CIRCULAR ORBIT 


PERTURB THE INITIAL CONDITIONS AT THE EARTH LAUNCH POINT BY 


ADDING A MANEUVER AT THE SPHERE-OF-INFLUENCE CROSSING OF 


THE EARTH AND MOON. 


THE COMPOSITE PATCHED CONIC WILL USE A TRAJECTORY 


THAT CAN BE ADJUSTED TO INTERSECT THE SPHERICAL 


MOON AT ANY DESIRED SURFACE IMPACT POINT ON THE 
NEAR-SIDE HEMISPHERE OF THE MOON. 





Figure 1.  Translunar Trajectory Design Flowchart. 
Precision translunar/transearth trajectory design requires time-consuming, step-by- 
step numerical solution of the exact differential equations of motion, including lunar and 


solar ephemerides. For preliminary design, the computation time can be reduced a 


hundred-fold or more by employing a closed-form "patched-conic" approximation to the 
precision trajectory. Since the patched-conic is a composite of the two conics, it can be 
formulated explicitly in terms of clementary functions of the key parameters. This 
provides better insight into design trade-offs than numerical computation. The patched 
conic approximates the gravitational effects of the Earth-Moon system. It ignores the 
small perturbation due to the Sun. 

Two approaches were considered for developing the patched-conic model. Both 
approaches employ Earth-centered and Moon-centered coordinate frames with respect 
to which the conic trajectories are formulated. Although the conic trajectories are planar 
with respect to thc Earth-centered and Moon-centered frames, they form three- 
dimensional curves as viewed in an incrtial barycentric frame. In the first approach, the 
coordinate frames are non-incrtial, being attached to the Earth and Moon which travel 
in (idealized) circular orbits (with respect to the Earth-Moon inertial frame). The second 
approach employs non-rotating, inertial Earth-centered and Moon-centered coordinate 
frames. These frames are attached, respectively, to a fictitious Earth and fictitious Moon 
which travel at constant velocity with respect to the barycentric frame. The implications 
of the second approach are briefly stated in Appendix A. The advantage of the first 


approach, pursued in this thesis research, is described in section A. (Wadsworth, 1991) 


А. TRANSLUNAR TRAJECTORY DESIGN 
The first consideration given to designing a translunar trajectory 15 to decide on 


what translunar trajectory types will be used. Since this thesis concentrates on fast Earth- 


return trajectories, eccentricities greater than one will be utilized. The reason for using 


non-inertial coordinate frames 15 discussed in section A.2. 


1. ‘Transhunar Trajectory Types 
Figure 2 illustrates all possible trajectory types that can be considered for 
Earth-Moon trajectory design. The radius of the Earth (R,), radius of the Moon (Ry), 
spacecraft true anomaly at launch (f), and spacecraft true anomaly at moon intercept 
(f) are projected in Figure 2. 
Since only trajectories with eccentricitics greater than 1.0 are being 
considered, only type I and II] hyperbolic trajectories will be used in thesis research (see 


Table 1). 


2, Non-Inertial Coordinate Frames 

The Earth-centered and Moon-centered coordinate frames are non-inertial 
frames, since they are in a state of gravitational free-fall about the Earth-Moon 
barycenter (center of gravity). As a consequence, the composite patched-conic provides 
a better approximation to the restricted, three-body spacecraft motion than would be 
suggested by simply joining two conics. 

Within the Moon’s region-of-influence, the terrestrial perturbing acceleration 
of the spacecraft is nearly the same as the centripetal acceleration of the Moon-centered 
coordinate frame. In this frame, these accelerations nearly cancel; this can be deduced 
from the facts that the distances of the spacecraft to the Earth and the Moon to the 
barycenter are comparable and the mass of the Earth is comparable to the equivalent 


mass at the barycenter. For the Earth-centered phase, there is a similar benefit over 


descending node of lunar orbit”. 


launch with latitude 


circle 
lunar orbit 


(direction of "Тә | 4 with direction 
earth rotation Е. / of motion 


indicated) 


ascending node of lunar orbit 
(a = 0) 





Figure 2. Equatorial plane projection. 


TRANSLUNAR TRAJECTORY TYPES: 


TYPE I posigrade : arc ab 

TYPE I retrograde : arc dc 

TYPE II posigrade : arc abc (e < 1 ) 
TYPE II retrograde : arc dcb (e < 1 ) 
TYPE III posigrade : arc dab 

TYPE II] retrograde : arc adc 

TYPE IV posigrade : arc dabe (e < 1) 
TYPE IV retograde : arc adcb (e < 1) 


Table I. Trajectory Types. 


most of the trajectory, especially for parabolic and hyperbolic trajectories. This 
cancellation does not occur in the case of patched-conics based on inertial frames as 
described in Appendix A (Wadsworth, 1965). 

The Earth-centered translunar conic trajectory may be cither elliptic or 
hyperbolic. The latter corresponds to shorter flight times and 1$ the only case that needs 
consideration for fast trajectories. For the same reason, the moon-centered conic 
approach trajectory is also modeled as hyperbolic in this study. Both the Earth-centered 
and Moon-centered frames are non-rotating. The motions of the Earth and Moon about 
their common barycenter are approximated by coplanar circular orbits. In general, the 
plane of these circular orbits does not coincide with the three dimensional trajectory of 


the spacecraft. (Wadsworth, 1991) 


В. SOLUTION FOR POINT-MOON INTERCEPT 
The least complicated form of a first order solution is for a point-moon intercept. 
Figure 3 illustrates the relevant parameters of the problem. X Appendix B, utilizing 
MathCAD 3.0, numerically solves the first order solution of point-Moon and massless 
spherical Moon intercept. In Figure 3 the following assumptions are made: 
e The kEarth-centered cartesian coordinate frame is non-rotating (indicated by 
A rand Z). 
* The Moon-centered cartesian coordinate frame is non-rotating (indicated by x,y,and 


z with the x-axis parallel to the X-axis, and the z-axis tilted at angle 1, for the Z- 
axis.) 


The initial point-Moon intercept solution represents a simple three-dimensional 
model; three-dimensional considerations are incorporated in the first-order solution solved 


in Appendix B. These considerations will be discussed in section B.2. 


p Three Dimensional Point-Moon Model 
Figure 3 illustrates the three-dimensional modcl of the first-order solution at 
point-moon intercept while Figure 4 further illustrates aspects of the three-dimensional 
model. Figure 5 delineates a flowchart of tasks to be accomplished in the three- 
dimensional model of the first-order solution of point-Moon and massless spherical Moon 
intercepts. This flowchart summarizes steps used to construct Appendix B. Symbols and 
abbreviations are defined in the Table of Symbols and Abbreviations. The hyperbolic 
trajectory yielded an energy requirement of 63.2 km"/sec? while the Apollo mission 
required a minimum energy requirement of 61.8 km"/sec. This results in a 226 increase 
in energy by using the hyperbolic trajectory, and as a consequence, will increase launch 
costs. Considering that the time of flight is reduced by more than 50%, the 2% energy 
increase penalty seems feasible for this mission. 
2. | Three Dimensional Considerations 
Figure 6 illustrates the relationship of orbit eccentricity and spacecraft true 
anomaly at moon intercept. Figure 7 is an exploded view of Figure 6 for user 
interpretation which shows eccentricities between 1.0 and 2.0. These diagrams must be 
referenced by the user to insure appropriate inputs for solution in Appendix B. Figure 
8 illustrates a cone sweeping out a circular area of the Earth. The Earth surface within 


the cone defines a zone of inaccessibility for Earth launch sites, given the specified true 


Earth Surface Launch Point 


for Reverse Trajectory Ascending Node of 


Moon's Orbit 





Figure 3. Three dimensional point moon model. 





Figure 4. Three dimensional model focusing on the point 
Moon. 
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Figure 5. Translunar Trajectory Program Flowchart. 


anomaly at moon intercept. This fact constrains inputs to the program in Appendix B. 


If the launch site is found to be in the restricted cone region a different trajectory must 


be used. 
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Figure 6. Eccentricity versus maximum and minimum true 
anomaly at Moon intercept. 
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Figure 7. Exploded view of Figure 6. 
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Figure 8. Three dimensional considerations. 


C. SOLUTION FOR MASSLESS SPHERICAL MOON INTERCEPT 

The final solution in Appendix B, is the massless spherical Moon intercept. An 
iterative technique is employed to obtain an impact on the lunar surface. This technique 
is illustrated and described in section C.1. Selenographic coordinates are utilized to 


obtain the solution for the massless spherical moon intercept and are described in section 


CZ. 


1: Iterative Technique for Lunar Surface Impact 
Figure 9 illustrates the mcthod to ascertain thc lunar surface impact point 


iteratively. The itcrative proccss is defined in the following steps: 
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* Determine point-moon intercept coordinates (R=pry, fy, Gy, and ty are inputs). 
¢ Select trial value of radial distance : Ry * pey-D. 


e Calculate the corresponding flight time from perigee, ty, and separation flight time 
В 


ST nd he nex iterauon: R, = Ret (Di - D)cosi,. 


> R,-D,+R,D, = cosa, and R, = R 
quantities). 


+(D = D)cosA„ı (bold print represents vector 


n-1 


* Continue until n 7 10 or | D,;-D|« 1 km. 


This iterative technique is employed in Appendix B, however, solutions are calculated 


manually by the program user without a terminal loop. 





P EM 
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SPACECRAFT 
TRAJECTORY 


TRACE OF P d 
SPHERE-OF-INFLUENCE 
OR LUNAR SURFACE 

( RADIUS D) 


LUNAR ORBIT 
(PROJECTION) 





Figure 9. Iterative technique model. 


2; Selenographic Coordinates 
Impact points on the near-side of the moon are expressed in selenographic 
coordinates, for which the zero degree longitude and latitude points at the Earth. The 
North, South, East, and West directions indicated on Figure 13 for latitude and longitude, 
are for an observer on the Earth. The lunar librations (Roy, 1988) are neglected in this 


definition. 


a. Euler angles 
Figures 10 and 11 express Euler angle transformations which are 
incorporated in Appendix B to relate selenographic coordinates to Earth-centered 
coordinate systems. The first Euler angle,(@), 1s zero because this angle is not required 
in this analysis. The second Euler angle rotation,(@), is the angle i in Appendix B. The 


third Euler angle rotation, ( V,,), is the angle $,, in Appendix B. (Goldstein, 1951) 


b. Lunar Impact 

Figure 12 illustrates a selenographic coordinate solution for lunar impact 
of the spacecraft. Since time reversal is employed in. Appendix B, the Moon moves 
backwards in its orbit for the translunar trajectory. The sign on o, is reversed in 
Appendix B for this reason. 

The selenographic latitude for surface impact is designated 6 
and the selenographic longitude for surface impact is designated œ. Figure 13 illustrates 
these latitudes, longitudes and their signs. The solution for æ has a factor of m added to 
it to insure that the positive X-axis points to the Earth. Figure 14 displays a locus of 


impact points for varying time offsets (Е). The axes units are in radians. Varying tetr 


14 


INCLLNATION OF LUNAR 
ORBIT PLANE WITH RESPECT 
TO EQUATORIAL PLANE 


SPACECRAFT 


EULER ANGLES : 
0 


1м 
Фм 





Figure 10. Second Euler angle rotation. 
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Figure 11. Third Euler angle rotation. 
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Figure 12. Lunar centered frame. 
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Figure 13.  Selenographic coordinates. 
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is equivalent to varying the input value of the lunar sweep angle (Įm). By providing an 
offset from the point-Moon intercept declination (6, и), Ше изег can target from North 


to South (see Figure 13). This factor can be combined with (tor) to produce an arbitrary 


location of the impact point on the lunar near side. 


(0$. ных 


$ 
750 seconds 0 


Figure 14. Locus of lunar impact points. 





D. PATCHED CONIC APPROXIMATION 
The method of calculation for the massless spherical moon intercept model can also 
be used for sphere-of-influence calculations if the radius of the sphere-of-influence is 


substituted for the radius of the Moon. The patched conic lunar impact point, using a 


finite lunar mass model, can then be approximated by an iteration process identical to the 
process used in Appendix B (Bate et al., 1971). 

From the final patched conic iteration, the impact point for the Moon-centered 
conic trajectory can be calculated. This calculation is made using the same method as 


Appendix B. 


15 


ПІ. APPLICATIONS 


Fast Earth-return trajectories have numerous applications, such as for requirements 
In evacuating a lunar base in an emergency or the establishment of a lunar-based 
Strategic deterrent force. This chapter will discuss such applications, their development, 


and technical feasibility. 


А. STRATEGIC PEACE INITIATIVE CONCEPT 
The ABM treaty prohibits nuclear weapons in space. If we look at a future world, 
this treaty might possibly be altered or not exist. The Strategic Peace Initiative is a 
hypothetical plan to place nuclear weapons on the Moon and to eliminate them from the 
Earth. The reasoning for this replacement is to prolong the strike capability and reaction 
time that now exists, thereby increasing stability. 
і. Background Information 
This section describes a hypothetical Strategic Peace Initiative, a plan to 
eliminate the risk of global nuclear war and simultaneously promote peaceful uses of 
space. The SPI concept could be a step toward reducing present nuclear confrontation 
as well as "denuclearizing" the Earth. (Wadsworth, 1958) 
Rather than using our current posture of SDI, the SPI concept can encourage 
peaceful exploitation of space (Gray, 1985). SPI resources could provide the logistics 
support for a new era of nonmilitary space activities, such as manned space exploration, 


servicing geostationary communication satellites, safe disposal of radioactive waste, clean- 
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up of hazardous space-debris, and eventually, deflection of Earth-threatening asteroids 
or comets. 

Although the basing of deterrent strike forces in space, particularly on the 
Moon, has already been proposed; this proposal is quite unique. The novelty of an SPI 
concept is that by a new treaty, the United States and Soviet Union would share the 
"high ground" of space. As with both the INF and Start treaties, success of the SPI treaty 
would depend on whether it serves the mutual interests of parties involved. (Wadsworth, 
1988) 

The free world’s offensive deterrent policies have prevented a global nuclear 
confrontation to date. Unfortunately, this policy is dangerously unstable because of the 
critically short decision time for launching a counterstrike to a perceived nuclear attack. 
SPI would eliminate or greatly reduce the danger of accidental nuclear war by increasing 
decision time from minutes, in the case of an SLBM attack, to a more comfortable two 
days. 

The primary justification for SPI is survival in a nuclear age, consequently the 
emphasis in this proposal is on deterrence. A secondary, economically-compelling 
justification is that much of the SPI investment in global security could also serve the 


civilian space program. They both require similar logistics elements: 


* A heavy-lift launch vehicle 
e A space tug (space/lunar transfer vehicle) 
* Earth/Moon space communications and tracking networks 


+ A manned lunar base 


Z: SPI Architecture 

The minimum-energy transearth strike trajectory has a 4.7 day flight time. 
Doubling the injection energy yields a three day trajectory. High-energy strike trajectories, 
taking two days or less, still permit nuclear strike weapons to escape the Earth-Moon 
system if aborted. Appendix B provides such a time of flight (approximately 1.9 days). 
Shortening the flight time from 4.7 days to 2 days results in only a two percent increase 
in transearth trajectory energy, which is equivalent to providing the additional velocity 
increment typical of an ICBM. Atmospheric re-entry velocity does not increase 
significantly between 4.7 day and 2 day trajectories, which simplifies re-entry vehicle 
design. To minimize vulnerability to intercept, the trajectories must be of the direct- 
ascent type without lunar parking orbits, in contrast to the more economical supply 
shuttle trajectories which must also be utilized. (Wadsworth, 1988) 

The SPI communications networks and surveillance systems would preferably 
be located at the natural svnodic satellite points, designated L1 through L5, which are 
stationary with respect to the Earth-Moon system and would minimize station-keeping 
fuel. These points are therefore an ideal location for placing these communications 
satellite systems. (Farquar, 1970) 

The two opposing lunar base zones, each approximately one hundred 
kilometers in diameter would contain the strike force missile silos and manned 
command/support complexes. These zones would be located on approximately opposite 


sides of the moon. They might occupy polar sites to take advantage of full-time sunlight 


for solar power. Determining optimum zone locations requires analysis if the optimum 
strike trajectory missions for both bases. (Wadsworth, 1988) 
Similar to the deployment of the stages of SDI, SPI would be implemented 
In stages and be subject to treaty compliance verification (possibly a United Nations role) 
at each stage. Both superpowers would deploy the following five stages: 
« NON-NUCLEAR DEPLOYMENT : Lunar bases, communications, and 
surveillance systems deploved. 


* FLIGHT TESTS: Test strike force with dummy warhcads to demonstratc accuracy 
and reliabilitv. 


* NUCLEAR DEPLOYMENT : Equip lunar bases with nuclear missiles and phase 
out Earth-based strategic nuclear forces. 


* VERIFICATION : Surveillance of lunar bases and Earth-Moon space; inspection 
and monitoring of Earth-based nuclear device manufacturing and test facilities; and 
inspection of translunar cargos. 

* SIRATEGIC REDUCTION EVOLUTION : Bilateral reduction of strategic 


nuclear forces to minimum required for astcroid deflections; incrcase in cooperative 
nonmilitary space ventures by utilization of lunar bases. 


В. SPI TECHNICAL FEASIBILITY 

Initial deployment of SPI relics on technology that has or will soon be proven by 
other programs. In fact, SPI can virtually be used as a spin-off of developed SDI 
technology. Programs such as TDRSS, NASA space station, OMV, NASP, HEDI, and 
the APOLLO lunar programs are examples of some of these programs. 

Challenges in technology do exist. These challenges include a precision guided re- 


entry vehicle, developing "smart" surveillance techniques, and prevention of adverse 
- £ 
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physiological effects of long-term residency in low gravity. At an interval of deployment 
of ten to twenty years, the fully deployed SPI deterrent force system probably would cost 
comparably to the proposed SDI system. (Wadsworth, 1988) 

Under the SPI treaty, a sovereign lunar base zone would be assigned to each 
superpower for basing its nuclear strike force. These zones would be located on 
approximately opposite sides of the moon. The optimum size and location of the zones 


would be a compromise among the following six basic requirements: 


e Minimize Surveillance Cost: Requires minimizing zone size; locate bases on the 
Moon’s visible side to take advantage of surveillance redundancy provided by Earth- 
based sensors to back up the spaced-based sensors. 


e Maximize Nuclear Survivability: Requires dispersal of hardened missile silos within 
a circular area of at least one hundred kilometers diameter. 


* Minimize Base Vulnerability: Requires maximizing separation between bases to 
deter pre-emptive strikes and reduce vulnerability of transearth strike forces to 
interception. 


* Minimize Transportation. Costs: Requires selection of base zone locations to 
minimize fuel cost for both supply shuttles which utilize parking orbits and strike 
forces which require direct-ascent trajectories. (This issue will be difficult to resolve, 
because if both strike lorces were to take advantage of the Moon's orbital velocity, 
then locations near the lunar poles will be favored. This induces a strategic 
problem, because the energy requirements for striking Northern-hemisphere Earth 
targets will differ between the polar bases.). 


* Avoid Base Overflight: Requires designing supply shuttle parking orbits not to 
overfly opposing base. 


* Minimize Communication Costs : Requires designing base communications on the 
near-side of the Moon. Far-side Moon bases require more communication system 
clements because its orientation with the Earth does not meet full-time line-of-sight 
communications link requirements. 


С. SPACE SUPPORT ACTIVITIES 

A large varicty of space support programs can be served by SPI. For example, 
lunar-based astronomy would be free of restriction imposed by Earth's atmosphere. A 
lunar base could be used as a staging platform for a manned expedition to Mars, taking 
advantage of a smaller round trip propulsion requirement for liftoff from the Moon, 
rather than from the Earth. The servicing of geostationary satellites could be 
accomplished with SPI resources. Transport of concentrated nuclear waste to a safe 
dumping site in the Moon would be unhampered by groundwater problems as on the 
Earth (Rosen, 1981). SPI resources could be used for debris removal from low-Earth 
orbits. This is an increasingly serious hazard to space missions. Also, safe disposal of 
space junk or dead satellites in higher orbits could be an SPI priority. 

Using the SPI surveillance network, the nuclear arsenal of the Moon can be used 
to deflect Earth threatening asteroids or comets. Collision with a ten-kilometer diameter 
asteroid has the explosive force of ten million one-megaton nuclear bombs. With an 
explosion of such magnitude, the impact dust cloud will plunge the world into a deep 
freeze similar to “nuclear winter” but much more devastating. Fossil records indicate that 
such an asteroid could account for extinction of the dinosaurs and cause for the Ice Age 


(Malove, 1985). 


D. SUMMARY 
The Strategic Peace Initiative Concept is a revolutionary plan for strategic nuclear 


deterrence. It is a plan which provides double value; a solution to accidental global 


nuclear war by increasing reaction ume to at least two days, and promotion of peaceful 
uses of space by developing permanently manned lunar bases and their infrastructures. 

The uniqueness of SPI is that it requires bilateral treaty agreement between the 
United States and the Soviet Union. Although the politics of such a treaty requires 
extensive thought, verification of the treaty can be accomplished. By looking to the 


future of peace, perhaps we can be at pcace with ourselves and with future generations. 
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IV. CONCLUSIONS AND RECOMMENDATIONS 


This final chapter summarizes the findings if this thesis. It also recommends 


research for further work concerning this topic. 


A. FINDINGS 
This thesis utilizes the selection of an efficient and operationally meaningful set of 
inputs and terminal parameters, to lead to a transearth or translunar trajectory targeting 


solution, which rapidly converges in an iteration process. 


В RECOMMENDATIONS 

This thesis lays the groundwork for an accurate patched conic method of designing 
fast transearth trajectories. Recommendations include the finalization of the patched 
conic approximation and the iterated impact point for the Moon-centered conic 
trajectory. A development of the SPI concept to include a cost/mission analysis would 


also be very useful. 
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APPENDIX A. PATCIIED-CONIC BASED ON INERTIAL FRAMES 


The patched conic trajectory can be formulated in terms of Earth-centered and 
Moon-centered coordinate frames which are inertial, being attached to a fictitious Earth 
and Moon which travel at constant velocity with respect to the barycentric frame. The 
Earth-centered conic trajectory approximation implies the lunar gravitational force is 
ignored so that the Earth is not accelerating with respect to the barycentric frame. 
Consequently, after the instant the spacecraft is launched, the center of gravitational 
force, the Earth, is replaced by a fictitious Earth which is imagined as moving in a straight 
line at constant velocity in the barycentric frame. In the meantime, the massless Moon 
continues on its circular orbit in the barycentric frame. In terms of the translating, Earth- 
centered frame, the lunar orbit appcars as a compound curve, rather than a circle. This 
fact must be taken into account іп targeting the conic trajectory for lunar intercept or 
crossing of the lunar-sphere-ol-influcnce. The same considerations apply to the Moon- 
centered conic phase where the Moon ts approximated by a fictitious Moon moving in a 
straight line at constant velocity in the barycentric frame. For greatest accuracy, the line 
is taken tangent to the Moon's circular orbit in the barycentric frame at the time of 


intercept. This requires an iterative solution. (Wadsworth, 1991) 


APPENDIENIS 


INPUT PARAMETERS: 


S/C dcclinalizs al eath launch point 


Gapul in deprecs) 


true cnomely of s/e ol more inlereest 


inpol in degrees) 


inclinalion of tie mosn orbit plone 


(inpol in degrocs : 


cecenlriody of e Uznsterer Urajectery 


\ 


(eccen >= еп } 
mn 


geocenlne allilude cf launch point 


(allude in kilometer) 


LUNAR TRAJECTORY PROGRAM 


Бов 
й 180 
бұ - 0.493928178314395 
radians 
Т . 


на 
M ' 180 


[M = 2.8623399732707 


radians 


д 


ham o 
M ' 180 


ing = 0.349065850398866 


redians 


eccenes 1.01 


26 


lonar sweep anele al mocn inleroepl M E 
180 
(inpul in degrees) 
Фм = 2.96705972839036 


radians 


poserade or relrograde ortil delerminalion 


poserode - +] 


relrogrede = -1 
«E opr 7! 
бм ОЙ :=0 ( nol iterated in program } 


INTERMEDIATE OUTPUT: 


ту 


бм = asin| sin. IM [sin] ФМ || 
д м = 0.059426145347819 


(spacecraft declination el moon ielercept) 





ам| = -0.986549223316463 


zl 


Q = 
M2 | cos! 6 M 


а м2 = 0.163464460888851 


самі | 
OW}? 





ам := ata 


ам = 2.97739131857372 


INPUT CONSTANTS: 


рем 7 384400 km : eath-moon mean distance 
Ге := 6378 km : earth mean radias 
гм := 1738 km : moon mean radins 


Ор = 7.292-107 rad/see: earth angular mlation male 
QM c= 152510 +— 
м 180 


о м = 2.661627109291352-10 ° rad/xc : 


moon meen orbilal rolalional rate 


2-1 
КЕ := 398600 km /see : earth gravilational constant 
4 9 
км := 4903 km /sec : moon gravitational constant 


OUTPUT VARIABLES end FORMULAS: 


ом = ом (dee lo the moon moving backwards in ils 


orbil for Lhe trandunar Urajeclory) 


ВЕ! TE thE 


Кє = 6.379*10° km : (s/c radial disdance al earth laench) 
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RMI ^ PEM 


RM1 = 3.844 10° km : (s/c radial distance al moon intercept) 


Км! 
КЕ! 


[ix 


(radial distance ratio) 


г = 60.26022887599937 


а r| 1eeccen*cos| fy ]]7 1 
eccen 


[Е = 0.726288833059664 


fem = Ем Е 
fem = 2.136051140211037 


бірл гЕМ `= sin{ TEM | NOTE: for S fem 


signppM = 0.84445214187218 if sen ts postive lher us +1 


И ten із negalive use -] 


Sfem ‘= | 


1 | ces[ à g | соз| дм || 


акмі 2 70.641452764362017 


Eu 


= a 25 Ў 
@EM2 7 J1^*EMI "Sfem** E OBT 


армэ = 0.767162532383019 NOTE : 
vsu (0емі | Ї | армі Т 
EM2 








is т'аЕ ОВТ 
if “ЕМІ > +} ther 


QEM = 2.267186781944015 


QE = AM FEM 

ap = 0.710204536629705 | ( 5/с пабі ascensa al carth Јевосћ ) 
бер г etl OM 

бЕ| = -0.338678825666912 

бЕ2= X fy 7 ÓM 

бЕ = 0.219826534971273 

(ы MAXI = зіп дє |і ём |+ со сем ]*cos[ 9g |*cos| à | 
[М МАХ! = -0.53563101113312 

ГМ_МАХ2 = 1-tM MAXI! | 


ЕМ_МАХ2 = 0.84445214187218 


ЧЕМ 


& EM 


M_MAX2 


(М MAX = 2 "аќап 





ам на 


ЇМ МАХ = 2.136051140211037 


И: Оор оро 


then Ім ЇМ MAX is constrained 
olherwise ase ар 


es calculated 


TE е аә әрәм) та 





| 


Aa = 1.570796326794897 – 0.496321 2320603151 
Cond rein: 


i Y 
aE MIN = ам 5 Аа 


aE yn = —0.164201335016073 » 0.496321232060315і 


ар med be greater than аЕмім 


п 

== 

| 
— 
— | 
аа) 


ессеп тіп” 


eccen min ” 0.97 4647 1757 23671 


n 
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| Skaris| 1t eccenss: ыл а: 
energy min ES ( = minimam specific energy 


al faenck } 


encrgy mig 7 -1.023795999276086 


ВЕ! = May ae ву 


cos] òg |*sin| fey | 





ВЕ! = 0.421410460300705 


sin| РЕМ | 





ВЕ? -| 


B £2 = 0.906870014913464 


Bp '- 2*atan 





1-8 J flight path azimath at levach: 


ВЕ? 
к>В>-к 


ВЕ = 1.135796261808675 


11 coe| à p |*sin| Bg | 


11 = 0.798478510707822 


12 ‚= MIU 


12 = 0.602023311789352 


ігі) (Lan donar trajectory inclination: 
04-1<4-я) 


і - 0.646032654504068 


aE сој = 0.699990269559789 


ап бє | 
tan(1) 





aE q `= 
aE Q? = 0.714152380463452 


1 ао] 


ао - зала 
aE o2 





ЧЕД = 0.795412455434794 


О := ар -аЕо ((гава опаг |гајескогу longitede- 
of- ascendinp- node) 


Q = -0.085207918805089 


СЕ! `= cos| дЕ | сов ав о, | 


СЕ! = 0.616325580024169 


35 





СЕ? := 





СЕ» з 0.787491447 196649 


1-Cry 





C p :- 2*atan 








СЕ? 


(sweep angle al leanch: 


СЕ = 0.906728172342672 УС у= 0) 
E 
| В кі Ї 1 +ессеп “05| Гр || 
криз joe (perigee distance) 


Eo 5.57010652841956 "107 


m kr '(1*eccen) 
PU | à 
| Кр (perigee velocity) 


Ур = 11.99319583652338 


2'kg (1- eccen)'kg 
Мене н а 
El p 


МЕ = 11.21107362995628 


(s/c velocity at lanach) 
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 Rp'(1-eccen) — (apogee disanco; valid only for ccce < 1 


B | - cccen 


ignore for eccen > = 1) 


В, = 5.57010652841956-10° 


(hyperbolic asymptotic Ure anomaly 
[asy 77 Eun - | 

balf- angle ; volid osly ог ессер 1 
[asy = 3.000756780023376 ignore for eccen < = 1) 
D :z 2*f asy 7x (deflection argie) 


D = 2.859920906456958 





: (direclrix-lo- focus distance ог axial offsel 
mu 
: of asymplolie from locus ; valid only for 


есееп > 1 ignore for ессев < = 1) 


4 
8 аѕу = 1.1085063487 25081 10 


REVO 


tr КЕ! “УЕ 


УЕ! = 0.934111286367402 


eccen*y p, зіп ГЕ | 


YE2 — 1 + eccen*cos| [E | 


7 £2 = 0.356981938872315 


УЕ := 2- мал (Night path angle al launch; 








-K 


a 
BU 


2 


ҮЕ = 0.365034946775507 


2°k С ессен)“ | | 
Ум := ыр ы кас (s/c velocity el moon intercept) 
вм; Вр 


Ум = 1.670175762283258 





ум = 9.1040525841947 18 


fMi = cos| ty | 


(ма = sin| fy ] 


, ессеп*у м1 “м2 


ум2 7 14 ессеп'ї мі ) 


у мо = 0.99457179716821 


1^ YMI 
YM2 


(fight path angle al moon ialercept) 





YM ‘= 2*atan 








ум = 1.466555060009538 


Вм! = 12: св ам -0)) 


Вм! = —0.600145987662734 











Вм> = 0.799890488437 212 


1 = 
Вм = зала a | (Might path azimuth el moon intercept) 
2 





BM 7 2-214479932657055 


Хм: Емі *cos| бм У сові ам | (RECTANGULAR COORDINATES OF 


THE NOON AF POINT-MOON 
Хм = -3.785601002578928 *10° 


INTERCEPT) 
Ум = Емј "сов] бм |'зпјам | 
Ум = 6.272482025241596 "10" 
2м := Вм: "3 6 м | 
7м = 2.282996752157727 *10' 
ЕЕ (Urea doner trajectory argument- 


w = 0.180439339283008 of- perigee) 


См := {м +0 


См 7 3.042779312553709 (s«cep angle al moon intercept) 











Хм: := e ian 7sin| à |1427<ф5| См-ум Ц'єсз| ам -Q]] ies 
+ SEAT 'sin ay - Q |] 
YM = s Era cin P и см ум oT j 





+ [1177м! "сої ам -0]| 


У м' = 0.048513058379955 


2 м' := Ум *12'с05| См ум | 


2. м' = -0.00545766973131 


[XM Xm || Ум'Ум' [2м 2м: | 


Км' = 
Км! 


(radial velocRy at poinl—-moon intercept) 


ВМ, = 1.648023680688986 
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15 
2*'r'Rp т | 
Пя (circular orbit period for radias R ) 


КЕ р 


D 4.137199589247244:10: 


| 1- cos[ f | 


2578 sin ГЕ | 


[Е| = 0.37999703209726 


1-соз| Гм | 


[M4 7 7.115369722384208 





Ep; ‘= Ш 





TEL - lee 
1 + eccen El 


E pi = 0.05361869041727 





Е р2 := 





Ер = 1.103630544045282 


4] 





1 1 
ессеп ~ 1 eccen - 1 
Е sx ros E “1 = 
shi | | + ессеп | | E 1 + ессеп а | | 


Еј = 9.05364438607 7804 
1 у -1 
ессеп - 1 eccen - 1 
E =i- | eÍ = | ET 
sh2 | | + ессеп МІ | | ' 1+сссеп М1 | 


Е аћ = 1-34171391602585 











_ [ эссеп*Е зы |-Ев1 
СЕТ с ар 
2*к*(ессеп- 1) 


tg, 7 0.089467283524084 


| eccen*E 2 |-E n2 


ма = 
2*к*(ессеп - ш 


tMi = 40.02754953820076 


E (eadh faunch fight lime relative 


to perigee passage) 


{р = 370.144008646907 
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м (poigl- moon intercept Might time 


relative lo pemgee passage] 


км = 1.656019615080179 -10° 


ЕМ = МЕ (poiet-moon lrandenar trajectory 


flight time) 
(ЕМ = 1.65231817499371 +10" 


INPUTS: 


Ко := 382986.2976625 (ary radiel disdance} 


OUTPUTS: 


2°КЕ| | (1- ессеп) "Кр | 


Мо = 
у Ro R 





р 


Уд = 1.6724659380983 


| Rp*(1 *eccen) | 


Ry:- Ro 


В. з 0.029233197611653 


2° atan 


== 
© 
її 





ај 


[eccen-1+R, | 








[0 = 2.861954051456415 
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Со = о +10 


C 9 = 3.042393390739423 
бо `= asin{ sin(i)*sin{ Cg ]] 


69 = 0.059657751937721 


ser, = Leu] 


Аа м1 = -0.996881642139325 


Аа м2 SM] 


tan(i) 





Aa M2 7 0.078911289214964 


Ла м := 2'аап 





Ла м = 3.062599237378809 


tan до | 


ee | sin(i) tan] C |] 


Aa | з -0.996857191524672 
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Гап] б0 | 
tan(1) 





Да » := 


Ac 2 = 0.079219566431173 





Аа ( = 3.062289992046826 


а, = ам“ дам 
а, г -0.085207918805089 


аб := Даџ+а, 


ао = 2.977082073241737 


_|Ер' Ур] 
| [ко "Уо | 
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¥ 1 = 0.104293660037493 


_ | eccen*y ; *sin[ fg || 


О | 1 + eccen*cos| fo | 


y 2 7 0.994546546158593 





cM 


ус -:- 2*atan 
ү 


Уд = 1.466312665340157 
By ‘= i2*cos| ag~ | 


B į = -0.600131267822715 





В) = 0.799901532303508 


a ама | 


2 





Bo = 2.214461530464985 


Хо := Ro'cos[ 59 |*ces[ a0] 


(RECTANGULAR COORDINATES OF 


-- — s 5 
Хо = -3.771433240671579 10 THE SPACECRAFT ON ITS TRAJECTORY) 


Хо В о"соз| бо | зіп/ ад) 


У 0 = 6.260990806241751 310" 


2.0 := Ко'з ёо | 


Zo = 2.283455104046146 *10' 
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ma: |1-сов| [0 |] 
sin| f | 


[о = 7.105421102050813 


|. (eccen - 1) 
а 


2 


pe оане 


| | (есеп = D 
(1*eccen) 


Е уз = 1.101755458701495 


“мо 
ме) 


1 I 
5; К (ессе 0) е И (eccen - 1) |. 
LE | Е] | EN | 


Е „з = 1.338578545032142 


(мо = 


(«мо = 39.8219787493874 


to “= Tp*tmo 


{0 = 1.64751474124978-10° 


[| eccen*E sh3 |- Ens | 


2*к*(сссеп- | 
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{өт = 9 (ime offset inpet) 


Arg ‘= tM ^to *tort 


Ar( 7 850.4873830398719 


Фм1 = cos| a M | соз| д м | 


фм = 70.984807753012208 


зіп à M | 
біп ім | 





9м2 “ 
Фм2 = 0.17364817766693 
Фм = 2 


$M 7 2.96705972839036 
$мо = мо М "Ато 


$ мо = 2.969323408665169 
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Х мо = РЕМ *cos| Фмо | (RECTANGULAR COORDINATES OF THE 
NOON AT PIKE t} 
0 


Х мо = -3.787102316833942 *10? 
У мо 7 PEM *cos| im |/зіл| Фмо | 
Y Mo = 6-191940094679928 -10° 


ZMO З РЕМ *sin| iy | зіп| ФМО ) 


2 мо = 2.253681886822944 -10“ 
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Do = ‚| [Хо-Хмо | *[ Уо-У мо | +[20-2мо ? | 
DISTANCE BETWEEN THE SPACECRAFT AND THE LUNAR CENTER AT TIKE | 
0 
D qo = 1-738000000264939°10° km 


| THE COAL I5 TO DRIVE THIS DISTANCE TO EQUAL THE RADIUS 
OF THE MOON FOR AN INPACT ON THE SURFACE OF THE MOON 


error ‘= D MOTIM 


error = 2.64939444605261 1°10 ^ 


[[Хмо - Хо |"Хо*[ Умо-Уо|"Уо*| 2мо- 20 |*2.0 | 


AES 
di |Ко'Рмо | 


А ро = 0.81263979792074 (= со [Хм] | 


R next = Ro+DMo*Apo 


( nomber to ase for next ileration of 
Ro } 


В лек = 3.843986656315015 10° 
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Х - Хм -Хмо 


Х = 150.1314255014295 


У: Ум-Умо 


У = 805.4193056166769 


Z:= Zy-Zyo 


Z = 293.148653347831 


( EULER ANGLE INPUTS ) 


д = ам 


ҰМ := м0 


WM 7 2.969323408665169 


5] 


OUTPUTS: 
(THE NATRIX STRUCTURE IS NUMBERED AS POLLOWS:} 


A, By Cy 
Ay By Cy ( ROTATION AND FRANSLATION } 
A3 B3 C3 


А у := совј У м |*cos(4) - cos(0) *sin( 9) *sin| Y y | 
А | = -0.985198313432347 

В| := -sin| Y y |“сов(Ф) - cos(0) *sin(9) *cos| V v ] 
В у = -0.171418444777856 

С} := sin(@)*sin(¢) 

C,=0 

A 2 `= соз ў м |*sin( ¢) * cos(0) "сов (ф)"зілі| Я м | 
А) = 0.161080647624348 

В) := 7sin| Y y |*sin(9) *cos(0)*cos($) *cos[ v y | 
В) = -0.925783585143099 

C 5 :7 -sin(0)*cos(4) 


C 2 = 70.342020143325669 
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A3 7 sin(0)'sin| Y M | 
Аз = 0.058628561051585 


Вз := sin(@) *cos| YM | 


Вз = -0.336957668364338 


С 3 := соз(д) 


C 4 = 0.939692620785908 


Aq B, C, 
до КА» во С» 
Аз Вз Сз 


-285.97295196007 
С =| -821.723449482487 
12.87940451383525 


Со = -285.9729519600736 


С, = -821.7234494824876 


Cy = 12.87940451383525 


X lun `= Co 
Y lun `= С, 
Z iun = С2 
Z 
5 := atan un 








2 2 
| Хуп + У уп | 


т SELENOGRAPHIC LATITUDE FOR SURFACE 


IMPACT ON THE NOON 


** 





s X lun 
a- =r ao --- пи 
2 2 
хь * Y lun | 


пн eee SELENOCRAPHIC LONGITUDE FOR SURFACE 


IMPACT ON THE NOON 
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